Non-adiabatic transport in a quantum dot turnstile 
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We present a theoretical study of the electronic transport through a many-level quantum dot 
driven by time-dependent signals applied at the contacts to the leads. If the barriers oscillate out 
of phase the system operates like a turnstile pump under a finite constant bias, as observed in the 
experiments of Kouwenhoven et al. [Phys. Rev. Lett. 67, 1626 (1991)]. The time-dependent cur- 
rents and their averages over succesive pumping periods are computed from the Keldysh formalism 
for tight-binding models. The calculation considers a sudden application of the pumping potentials 
at t = which leads to transient features of the time-dependent and averaged currents during the 
first pumping cycles which turn out to be important in the high-frequency regime. We show that 
in the transient regime the efficiency of the system as a pump is rather poor because it mainly 
absorbs charge from both leads in order to fill the levels located below the bias window. Under a 
finite bias and a low-frequency pumping signal the charge transferred across the system depends on 
the number of levels located within the bias window. The internal charge dynamics and the role of 
energy sidebands are investigated. The so called satellite peaks of the averaged current are observed 
also in the transient regime. 



PACS numbers: 73.23.Hk, 85.35.Ds, 85. 35. Be, 73.21. La 



I. INTRODUCTION 

The ability to control the transport properties of 
semiconductor quantum dots by time-dependent per- 
turbations (e.g. microwave signals or optical pulses) 
allows the observation of photon-assisted tunneling, 
charge pumpingLSiS and coherent Rabi oscillations. 4 
Also, pump-and-probe techniques were used to estimate 
relaxation rates and to control spins in quantum dots. 5 
A common point of these experiments is that the time- 
dependent driving potential is applied on the system it- 
self, i.e. on a central metalic gate defining the quantum 
dot. 

Some time ago Kouwenhoven et al& proposed a differ- 
ent setup, in which a quantum dot is coupled to source 
and drain reservoirs by oscillating tunneling barriers. 
Technically this is achieved by applying radio-frequency 
signals to the metalic gates that control the opening of 
the quantum dot to its surroundings. The two barriers 
at the contacts are varied in such a way that the system 
undergoes a cyclic transformation and, under a constant 
bias applied on the leads, an integer number of electrons 
is transmitted during one cycle. Therefore the system op- 
erates as a turnstile pump. An important feature of the 
turnstile configuration is that the pumped current has a 
definite direction due to the finite bias. We remind here 
that originally the concept of parametric charge pumping 
was introduced in the context of a net current genera- 
tion in an unbiased system^ As noted in the literature, a 
symmetry breaking is necessary in order to get a nonva- 
nishing current without a bias. In spite of the fact that 
the turnstile operation was experimentally observed some 
time ago, it attracted little attention in the theoretical 
literature (see Refs. Ilfjl3lll4lll7l below). The purpose of 
this work is to explore the transport properties of turn- 
stile quantum dots submitted to time-dependent signals 



of arbitrary amplitude and frequency. 

At theoretical level the quantum pumping was dis- 
cussed basically within two frameworks: the adiabatic or 
Floquet scattering theoryi^il and the non-equilibrium 
Green-Keldysh formalism i 14 ' 15 i 16 i 17 The scattering ap- 
proach was primarily designed to describe the adia- 
batic pumping, in which the driving potential varies very 
slowly. The timescale on which the applied signal varies 
significantly exceeds the time needed for the electron to 
pass through the system. The key result of the adiabatic 
scattering is a current formula in terms of an instanta- 
neous (frozen) S matrix. This matrix is computed per- 
turbatively and only the linear term in frequency is usu- 
ally retained. A rigorous mathematical treatment of adi- 
abatic quantum pumping^ recovered the BPT formula 
given inii Within this framework the relation between 
resonant transmission and quantized pumped charge in 
an unbiased turnstile was analyzed* 1 ^ An extended scat- 
tering formalism for studying both adiabatic and non- 
adiabatic quantum pumping was developed in Ref . and 
uses the Floquet theory and an S matrix depending on 
two-energies. Recently Mahmoodian et al. 13 have com- 
puted the stationary current for a quantum wire submit- 
ted to alternating 5-like voltages. It was shown there that 
the current displays multiphoton peaks as a function of 
the Fermi momentum of the leads. 

Using the Green-Keldysh formalism Q-f Sun and T. S. 
Linii have computed the current through a single level 
quantum dot when rectangular or harmonic potentials 
are applied at the contact to the leads. As it is well 
known, this model is exactly solvable within the wide- 
band limit (WBL) approximation since the Dyson equa- 
tion for the retarded Green function is greatly simplified 
due to the leads' self-energy which within WBL is simply 
a delta-function.— Later on Wang et al. used the Keldysh 
approach to investigate the non-adiabatic charge pump- 
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ing in the presence of photon-assisted tunneling. The 
explicit calculation was done within WBL and for an un- 
biased double barrier pump driven by a local sinusoidal 
signal. It was shown that at large frequencies a nonvan- 
ishing current is generated even with a single-parameter 
perturbation. Moreover, a sign change of the pumped 
current was reported when the Fermi level of the leads 
crosses the eigenvalues of the system. This feature was 
predicted also by Buttiker and Moskalets*S 

Further progress was achieved by L. Arrachea for tight- 
binding models and periodic potentials^ The method 
developed in this paper allows the calculation of the 
d.c. component of the pumped current once the par- 
tial Fourier transform (i.e. the Fourier transform with re- 
spect to one time only) of the Green functions is known. 
It was shown that the d.c. component can be related 
to a transmission function T(u>) which is interpreted as 
the difference between the probabilities of tunneling out 
from and into the system. The numerical simulations 
are performed for unbiased one-dimensional pumps and 
the pumping potential is described by a diagonal time- 
dependent term added to the energy at the contact sites. 
For two harmonic potentials the retarded Green func- 
tion is computed perturbatively for weak pumping am- 
plitude and pumping frequency. The connection between 
the Floquet scattering and the Green-Keldysh function 
formalism for time-dependent transport was discussed in 
Ref0. 

Besides the scattering theory and NEGF approach the 
pumping problem can be addresed via time integration 
of the Schrodinger equation by the Crank-Nickolson ap- 
proximation, as proposed by Stefanucci et al. in a series 
of papers i 19 i 20 Their setup starts from the ground state 
of the unperturbed but coupled system and has there- 
fore the advantage of introducing naturally the bias as 
a perturbation. We remind that the Keldysh formal- 
ism requires a partitioning of the system into 'central 
region' and 'leads', the perturbation being the coupling 
between them which is established usually adiabatically 
in the remote pastil The Keldysh approach is however 
appropriate for studying the transport in the turnstile 
configuration which requires to connect and disconnect 
periodically the pump from the leads. 

Recently the equation of motion method was applied 
to compute various currents in ID pumps coupled to fi- 
nite wires with constant chemical potentials* 2 ^ Such an 
assumption is questionable for finite systems, especially 
when one is interested in the long time behavior. 

In this work we are primarily interested in the tran- 
sient effects on the transport properties of a many-level 
quantum dot turnstile which is submitted to pumping po- 
tentials at t = 0. We believe these effects could seriously 
affect the transport properties of nanostructures driven 
by fast oscillating signals. Also, we mention that most of 
the previous theoretical approaches present calculations 
of the stationary averaged current which implies either 
to look at the long-time limit behavior—, either to con- 
sider small pumping frequencies ^ Since in the stationary 



regime the transient effects are presumably washed out, 
one does not need to specify how and when the pumping 
signal is turned on. Note also that for a proper appli- 
cation of the Keldysh formalism it is crucial to have a 
well-defined equilibrium state of the decoupled system. 

While in the adiabatic regime the main advantage is to 
express transport quantities in terms of the frozen scat- 
tering matrix up to errors of 0(u> 2 ) or even 0(o> 3 ) (see 
the higher order corrections in Ref. @), the NEGF formal- 
ism covers the entire frequency range, allowing therefore 
an equal footing treatment of adiabatic and nonadiabatic 
pumping. The Green functions needed in the current for- 
mula are computed using a recently developed method 23 
which solves the integral Dyson equation exactly by trans- 
forming it into an algebraic equation. Through this pro- 
cedure the Green functions are computed taking into ac- 
count all back-and-forth scattering processes. 

The content of the paper is divided as follows. Section 
II describes the model and gives the relevant equations, 
as well as the considered pumping potentials; more de- 
tails about the formalism should be traced back from Ref. 
[23l . Section III is the main part of the paper and presents 
the numerical results and their discussion. Conclusions 
are summarised in Section IV. 



II. THE MODEL 

Within the tight-binding model which is adopted 
throughout this work the Hamiltonian of the system con- 
tains three terms: the semiinfinite leads (Hl), the quan- 
tum dot turnstile (i?s) an d the time-dependent pumping 
signals ifr(^): 

H(t) = H s + H L + H T {t). (1) 
H$ has a usual tight-binding form 

N 

Hg (fm + V g )d\ n d m + 22 t mn dl n d n . (2) 

m=l (m,n) 

Here t mn are hopping terms, (m, n) denotes nearest- 
neighbor summation over the system sites. e m is the on- 
site energy and the diagonal term V g simulates a plunger 
gate potential applied on the system. N is the number 
of sites in the dot. 

Following the experimental setup from Ref. @ we de- 
scribe the oscillating tunneling barriers between the dot 
and the leads by time-dependent hopping terms (l,r de- 
note the left and the right lead): 

H T (t) = V a {t)(c\j mci + h.c). (3) 

a— l,r 

Here Ci a and c\ denote the annihilation/creation opera- 
tors on the z-th site of the lead a which is connected via 
the nearest neighbor hopping V a to the site m a of the 
dot. Similarly, d ma and d/ m correspond to the site m a 
of the dot which is coupled to the lead a. 
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FIG. 1: (Color online) The pumping potentials applied on the 
left (Vi - solid line) and right (V r - dashed line) contacts. We 
show two pumping cycles k = 2. The parameter A = 0.6, 
vi = v r = 0.75 and the frequency u = 1. 

The signals applied at the contacts between the dot 
and the leads have a trapezoidal form and are defined 
with the help of a function with period 2tt introduced as 
follows (A is a positive number): 

!1 if LOt £ [0,7T- A], 

l-£(wf-(7r-A)) if «*G[7r-A,7r], , , 

if LUt e [7T,27T-A], W 

±(ut- (n - A)) if cot g [2tt - A,2tt], , 



In the above formula the retarded and lesser Green func- 
tions are given as usual in terms of Heisenberg operators 
G R (t,t') = -i«(t-t0({c il (f),4 l (t)}) and Suit,*) = 
{(c^ (t')d nl[ (t)). fi(E) is the Fermi function of the left 
lead, tL is the hopping energy on leads and T l is the 
linewidth depending on energy and time: 

r l (E;t,s)=p(E)V l (t)V l (s), (8) 

containing the pumping potentials at different times and 



G R {t,t') 
G<(t,t>) 



where G R ' A (t,t') are the retarded and advanced Green 



Then the pumping potentials of period T are given by 
the relations: 

Vi(ut) = vif(LJt) (5) 

V r (ut) = W P (l-/(wt)), (6) 

where to = 2n/T is the frequency and vi jr are the am- 
plitudes of the pumping signals. It is useful to introduce 
the number of pumping cycles k considered in the nu- 
merical simulation. We show for clarity in Fig. I a train 
of two such pulses (i.e. k — 2) that we use to simulate 
the turnstile configuration. The quantum dot is coupled 
suddenly to the left lead at to = while the right con- 
tact is off. In the range [kT/2 — A, feT/2] the sample is 
simultaneously isolated from the left lead and connected 
to the right lead. This switching is done linearly at a 
slope 1/ A (note that a larger A implies a slower onset of 
the couplings). In the second halfperiod of each pumping 
cycle the right contact is open, allowing thus the charge 
pumping. The cycle is completed by lowering the left 
tunneling barrier (i.e. increasing Vi(t)) and turning off 
the coupling to the right lead. 

As it is widely known, the standard application of the 
Keldysh formalism leads to the following formula for the 
current entering the system from the left lead (here we 
take for simplicity one dimensional leads; a many-channel 
formula and more details are to be found in Rcf. 23): 



(7) 

I 

the density of states at the endpoint of the semiinfinite 
one-dimensional lead p{E): 



p(E) = 6{2t L -\E\)¥-!jp . (9) 

The retarded and the lesser Green functions are com- 
puted from the Dyson and Keldysh equations: 

(10) 
(11) 

functions of the isolated dot and E K,< are the retarded 



Ji(t) = ~lm( / dE / d S e-* E ^T l (E;t, S )(G? l (t, S )f l (E)+G<(t, S ))). 

n J -2t L JO 



= G$(t,t') + / dhG^h) / dt 2 -Z R (t u t 2 )Gn{t 2 ,t') 



t f t' 

dtiG R (t,t{) / dt 2 ^ < (t ll t 2 )G A {t 2 ,t l ), 
o Jo 
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and lesser self- energies. It worths mentioning that both 
self-energies contain the known Green functions of the 
semiinfinite leads but they are also quadratic functions 
of the pumping potentials whose time variable is differ- 
ent (see Eq. [8]) . Therefore, their time-dependence is 
much more complicated than in other approaches were 
the pumping signals are applied to the system or to the 
leads and are described by diagonal terms in the Hamilto- 
nian. In particular, the algorithm taken in Ref. would 
we difficult to use. The time-dependent occupation num- 
ber can be computed from the lesser Green function of 
the dot: 



JV 

N(t) = Im Gl m (t, t) = Y, N "^- ( 12 ) 



N 



N m (t) are on-site occupation numbers and will be used in 
the next section to gain information about the internal 
charge dynamics during the pumping cycle. A similar 
formula can be written down for the current J r (t) flowing 
from the system towards the right lead. 

As we have said, taking explicitely into account the 
initial instant when the pumping signals are applied leads 
to transient effects.— One consequence is that the period- 
averaged currents depend on the period index k. We 
introduce therefore a fc-indexed period-average for the 
currents (the fc-th period covers the interval \tk-i, tk] and 
to = 0): 



J a, 



i 

f 



dtJ a (t), 



a 



Lr. 



(13) 



Although the approach taken in this work does not in- 
clude the electron-electron interaction it captures the ba- 
sic known features of turnstile pumps: the quantized 
pumped charge in low frequency regime and the satellite 
peaks due to photon-assisted tunneling in high-frequency 
regime. Moreover, most of the results are presented at 
rather strong coupling to the leads when the dot is fairly 
open and the Coulomb blockade effects are not impor- 
tant. In a recent work Splettstoesser et at— proposed 
a method for dealing with Coulomb interactions in adi- 
abatic quantum pumps. This approach is based on the 
quantum Master equation and uses a perturbative expan- 
sion in the tunnel coupling. In the case of the unbiased 
quantum dot turnstile these authors found that one has 
to consider the second order term in the tunnel coupling 
to get a non vanishing pumping. 



III. NUMERICAL SIMULATIONS 

In this section we present the main numerical re- 
sults and discuss the transport properties of the turn- 
stile pump in different regimes. The bias, the frequency, 
the energy, the hopping constants on the leads, the cou- 
pling strengths and the gate potential will be expressed 
in terms of the hopping energy of the central region t o 
which is chosen as energy unit. 
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FIG. 2: (Color online) (a) The pumped current in the right 
lead (J r ) and the current entering the system from the left 
lead (Ji) for a 2-site turnstile, (b) The occupation number 
of the dot N(t) and the on-site occupations Ni(t), i — 1,2. 
The pumping potentials are also given. During the second 
halfperiod of each pumping cycle the turnstile expells one 
electron to the right lead. The bias is fixed to W = 3.0, 
io = 0.6, vi = v r = 0.75 and kT = 0.0001. 



The current is therefore given in units of eto/Ji and the 
time expressed in units of Ti/t p. We take e — h = 1. The 
bias window (BW) is defined as the difference between 
the chemical potentials of the leads W — /i; — fi r . To 
make a connection to physical units one could take for 
example the energy unit as tjj = 0.1 meV. Then the 
frequency unit would be lo ~ 25 GHz and the time unit 
t r~ 7 ps. The termal energy kT = 0.0001 in all numerical 
simulations. 

We consider first a two-site one-dimensional turnstile 
submitted to a finite bias and modulated by the trape- 
zoidal signal introduced in Section II. In order to simulate 
the conditions of the experiment performed by Kouwen- 
hoven et al& we set the bias window to W = 3.0 with 
respect to the chemical potential of the right lead /i r = 
such that the highest level of the isolated dot E\ = 1 is 
located in it. Fig. 2a shows the current Ji{t) from the left 
lead towards the turnstile and the current J r (t) pumped 
into the right lead during two pumping cycles. The fre- 
quency is lj = 0.6 and the maximum height of the tun- 
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neling barriers is V\ = v r = 0.75. As expected, a nonva- 
nishing current is generated in the right lead during the 
second halfperiod of each pumping cycle. The pumping 
mechanism is proved by the behavior of the occupation 
number N(t) which is given in Fig. 2b. At the begin- 
ning of the first cycle the system collects charge from 
the left lead and since the coupling to the right lead is 
zero there are almost two electrons in the system after 
one halfperiod. We shall call this halfperiod the charging 
halfperiod. The charge dynamics within the system is 
also given in Fig. 2b through the occupation numbers of 
the two sites Ni and N2 . For the clarity of discussion we 
have also included in the figure the two potentials applied 
on the leads. The first site is rapidly populated up to 0.5 
and then stays at this occupation while the second site 
starts to be filled. The current J; decreases in this short 
time range. 

Then we see that N2 increases faster than Ni and both 
are reaching a constant value. Due to the small frequency 
considered here the halfperiod of the pumping cycle ex- 
ceeds the time needed for the system to be completely 
filled with electrons and therefore, as long as the cou- 
pling to the right lead is still turned off, the total occu- 
pation number is steplike. The current J\ becomes very 
small in this range because as both sites are filled it is 
more difficult to inject charge (notice the lower slope at 
which Ni increases). The pumping starts effectively in 
the second half of the cycle and leads to the transmis- 
sion of one electron in the right lead, as observed in Ref. 
6. During the pumping halfperiod the occupation of the 
right contact drops quickly to 0.75 but then decreases at 
a lower rate as the first site is depopulating. Note that in 
the pumping sequence N\ drops more rapidly than iVa- 
The occupation number at the end of the pumping cycle 
does not decrease below 1 because the lowest level of the 
dot Ei = — 1 is well below the bias window and cannot 
contribute to transport. 

From the above discussion it is clear that the efficiency 
of the pump depends on two facts: first, the charging 
halfperiod should allow the complete filling of the turn- 
stile and second, the levels within the bias window must 
be entirely depopulated during the pumping halfperiod. 
In the intermediate regime to ~ Vi >r (not shown) the 
system still transfers one electron, the difference being 
that the occupation number decreases faster and more 
abruptly than in Fig. 2. The transient effects on the time- 
dependent currents and pumped charge that appear as 
the pumping frequency increases are captured in Figs. 3a- 
d. At uj — 1.25 the shape of the currents is similar to 
the one in Fig. 2a. However, the current pumped dur- 
ing the first cycle is smaller that the ones corresponding 
to the next cycles and the occupation number plotted in 
Fig. 3d shows that less charge is pumped. Also, there 
are no more steps in N(t). The behavior of the occupa- 
tion number helps us to identify the new features of the 
charge dynamics in the nonadiabatic regime. 

On the other hand, the coupling to the right lead closes 
before the level within the bias window depopulates. As 
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FIG. 3: (Color online) The shape of the currents J r and Ji 
depends on the pumping frequency: (a) to = 1.25, (b) u = 
1.75, (c) uj = 3.5. At large frequency J r takes negative values 
in the first pumping cycles, (d) The occupation number N(t) 
corresponding to the three frequencies taken in (a), (b), (c). 
W = 3.0, vi=v r = 0.75, kT = 0.0001. 
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FIG. 4: (Color online) (a) The average currents J Tt k and Ji t k 
as a function of the period index k. There are three pairs 
of curves, corresponding to the frequencies ui = 0.6 (k = 4), 
uj = 1.75 (fc = 5) and uj = 3.5 (k = 11). For each frequency 
the two averaged currents are drawn with the same type of 
line (color); J r (Ji) has a positive (negative) slope, (b) The 
average occupation number (see the discussion in the text). 
W = 3.0, vi = v r = 0.75, kT = 0.0001. 



a consequence, at the end of the pumping cycle there is a 
residual charge (~ 0.1) which is stuck within the bias win- 
dow. By further increasing the frequency up to ui = 1.75 
another effect appears in Fig. 3b. Although there is a 
nonvanishing pumped current during the first cycle J r 
takes negative values at the end of the cycle, which means 
that in this interval the system absorbs charge rather 
than pumps it. Looking at the corresponding occupation 
number in Fig. 3d one infers why the system does not act 
like a pump over the entire pumping sequence. N(t) goes 
slightly above 1 then drops to 0.9; during this interval the 
system pushes electrons to the right lead and therefore 
J r (t) > 0. Then the occupation number increases, lead- 
ing to negative values of J r (t) (the pumping period in 
Fig. 3b is 3.7). The physical picture behind this is the 
following: i) In the first halfperiod both levels are popu- 
lated, though not completely; ii) During the first part of 
the pumping sequence the level within the bias window 
depopulates, generating therefore a positive current J r ; 
iii) When this happens the total occupation number 



FIG. 5: (Color online) (a) The occupation number at differ- 
ent values of the bias window and different allignment of the 
turnstile levels with respect to the bias window, (b). The 
pumped current J r (t) corresponding to the occupation num- 
bers in (a). For all curves the frequency is set to u) = 0.3. We 
show only the first pumping cycle. 

goes below 1; also, since the chemical potential of the 
right lead is higher than the lower level of the dot, the 
system absorbs charge form the right lead and the occu- 
pation number increases again. 

Finally, Fig. 3c shows the currents in a highly non- 
adiabatic regime w = 3.5. The steplike features within 
each pumping cycle are washed out and on the first cycle 
the system does not pump any charge (the occupation 
number is simply increasing and stays below 1). More- 
over, even when a stable pumping regime is achieved (for 
t > 5) it is not effective at all since very little charge is 
pumped. 

The transient regime should be noticed as well in the 
period- averaged currents. In order to check this we give 
in Fig. 4a these currents for the three frequencies consid- 
ered in Fig. 3. The following things are observed: i) In 
all three cases the d.c. components of the currents be- 
come eventually equal and their value does not depend 
anymore on the period index. Moreover, the averaged 
current is conserved. The passage to this 'stationary' 
regime is faster at low frequencies; ii) In the transient 
regime J;^. exceeds J r ^ since there is a net charge ac- 
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cumulation in the dot along each transient cycle (see the 
occupation numbers in Fig. 3d); iii) At large frequencies 
J, % k takes negative values in the transient regime (see the 
first cycle at uj = 3.5) because the system absorbs charge 
from the right lead. 

Fig. 4b shows that the average occupation number 
depends strongly on the period index in the transient 
regime and settles down to 1.5 in the long-time limit. 
Notice that this means that the level which contributes 
to the transport is half occupied. 

The bias applied on the leads is an important param- 
eter in turnstile operation because it controls the num- 
ber of levels giving the main contribution to the cur- 
rent. As reported experimentally^ the number of elec- 
trons pumped during one pumping cycle is given by these 
levels. Fig. 5a shows the occupation numbers of a three 
sites turnstile pump submitted to a different bias. We 
take a small frequency (uj — 0.3) and therefore the pump- 
ing cycle is long enough to allow the filling of the low- 
est levels (the dot has three eigenstates E± — ±1.4 and 
Eq = 0). For /i; = 3 and /z r = 1 only the highest level is 
located within the bias window and the pump transfers 
one electron. The occupation number at the end of the 
cycle is N — 2. We emphasize that the frequency is much 
smaller than the gap between Eq and E_ and therefore 
the lowest level cannot give an important contribution to 
the current via excited sidebands. 

By decreasing the chemical potential of the right lead 
to jji r = the middle level alligns to fi r and the turn- 
stile pumps 1.5 electrons at W = 3. Note that due to 
the rather large coupling to the leads there is no charge 
quantization condition to be fullfilled. The dotted line 
shows the occupation number for the same bias W with 
Hi = 1.5 and fj,i = — 1.5. In this case we have a sym- 
metric bias window (as marked in the figure) contain- 
ing the middle level in the center. The other two levels 
are close to the chemical potentials of the leads. The 
transferred charge is Q p ~ 2, suggesting that each level 
alligned to one chemical potential pumps only half elec- 
tronic charge. This happens because the level close to 
fi r (pi) are difficult to depopulate (populate). Also, by 
tuning the width and the position of the bias window it 
is possible to transfer the same charge in different ways 
(i.e. by involving different levels of the pump). For ex- 
ample, seting fii — 1.5 and fj, r = one can still transfer 
one electron at a bias window W = 1.5. Of course, the 
efficiency of the pump increases when the number of lev- 
els participating in transport increases. Fig. 5b shows the 
pumped currents associated to the occupation numbers 
in Fig. 5a. Although both cases W — 1.5 and W — 2 
correspond to a pumped charge Q p — 1, the current is 
higher at W — 1.5. Note also that the currents have a 
peak structure when two levels participate in transport. 

The discussion around Figs. 3 and 4 suggests that at 
a given frequency the transient regime would cover more 
pumping cycles if the number of levels located below the 
bias window increases. This is confirmed in Fig. 6 which 
shows the occupation number of a 3-site turnstile. We 



take fii = 3 and /i r = 1 and in this case the BW contains 
only the highest level E+ = 1.4. 




FIG. 6: (Color online) (a) The total occupation number N(t) 
and the on-site occupation numbers Ni (i, 1, 2, 3) for a 3-site 
turnstile submitted to a bias W = 2 and to a pumping signal 
of frequency u — 1.57. 

N(t) shows clearly that the other two levels are filled 
only after k = 5 pumping cycles. Notably, even in the 
transient regime the system pumps a small amount of 
charge (~ 0.25) because of the highest level from the BW. 
This level cannot be completely filled in the transient cy- 
cles because most of the charge populates the lowest two 
levels. The behavior of the individual occupation num- 
bers Ni, i = 1,2,3 confirms the intuitive picture of the 
internal charge dynamics. It is easily seen that for any 
pumping cycle except the first, the occupation number of 
the 3rd site decreases even in the first halfperiod, namely 
during the charging process. Since in this time range the 
contact to the right lead is turned off the charge from 
the right contact can only flow back to the middle site. 
Indeed, N 2 increases as N 3 decreases. 

As the dot charges from the left lead the inverse flow is 
noticed. In the pumping halfperiod the middle site occu- 
pation firstly increases because the charge accumulated 
in the first site passes to the right contact which is now 
open, and finally decreases. All these internal bounc- 
ing trajectories are taken accurately into account in our 
calculation because the method we used to solve for the 
Dyson equation gives the entire matrix of the Green func- 
tion (not only the diagonal matrix elements entering the 
current formula). 

We turn now to another important feature of time- 
dependent transport, namely the contribution due to in- 
elastic scattering processes in which the incident elec- 
trons gain or loose energy quanta from the driving fields. 
Within the Floquet scattering approach it was shown&S^ 
that the outgoing electrons can have any energy E n = 
E + nhu> in the so called sideband ladder, E being the 
energy of the incident electrons. In the present model 
the pumping potential is not periodic because of the con- 
dition Vi(t) = V r (t) = for t < 0, therefore the Flo- 
quet theorem cannot be rigorously applied. At best, one 
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FIG. 7: (Color online) For a 3 site turnstile with all the levels 
above the bias window one can notice sideband contributions 
to the two time-dependent currents at frequency (a) u> — 1.0, 
(b) uj — 1.6 and (c) uj = 2.75. (d) The occupation number at 
different values of the frequency. 



can use results of the Floquet theory in the long time 
limit where all relevant quantities will oscillate with the 
frequency of the driving signals (see Ref. [H for a dis- 
cussion). Nevertheless, inelastic tunneling processes are 
physically possible especially in the transient regime and 
should be noticed in the averaged currents. In particular 
one expects to see additional contributions to the current 
at frequencies matching the gaps ± Ei between the 
levels of the isolated system and the chemical potentials 
of the leads. 

In order to check these features we tune the bias win- 
dow below the three levels of the 3-site pump by setting 
fil = —2.5 and \i r = —3.5 such that E- — Hi ~ 1. We 
take a small pumping amplitude vi = v r = 0.35 and 
look at the currents for several frequencies, as shown in 
Figs. 7a-d. For small frequencies (uj = 0.3 and uj = 0.5 
- not shown) there is no pumped current and J/ takes 
both positive and negative values, while N(t) increases 
and decreases accordingly during the charging halperiod. 
This suggests that the system repels incident electrons 
back to the left lead. This happens because on one hand 
there is no level below or within the bias window and, 
on the other hand, the frequency is too small to allow 
the population of the lowest level. The situation changes 
for uj = 1 (see Fig. 7a) . A pumped current appears in all 
three cycles presented. We related this to the occupa- 
tion number plotted in Fig. 7d which shows that there is 
a charge pumped into the right lead. When comparing 
Fig. 7a and Fig. 7d (the dotted line) one notes that in the 
second half of the charging cycle <// < and N(t) de- 
creases. This signals the energy relaxation process from 
the lowest level of the system to the left lead only, be- 
cause the contact to te right lead is not yet open. Only 
the second (lower slope) decrease of N(t) is associated 
to pumping trough the right lead. At us = 1.6 (Fig. 7b) 
the pumped charge and current increase. J; is mostly 
positive because the fast oscillating signals prevent relax- 
ation to the left lead. The pumped current still shows a 
peak structure. In Fig. 7c we show a highly non-adiabatic 
regime at uj = 2.75. In this case the second level Eq = 
will contribute as well. In order to avoid tunneling from 
the right lead to the level E_ we have considered for this 
curve fi r = —5. The peak structure of the pumped cur- 
rent dissapears and a periodic regime establishes slowly 
after 9 cycles. As we have said in the introduction, the 
present approach takes into account all tunneling pro- 
cesses between the leads and the system and within the 
system. This nonperturbative treatment might be crucial 
for describing the relaxation processes mentioned above, 
especially in the transient regime where the resonances 
are not well defined. 

Finally we investigate briefly the satellite peaks ap- 
pearing in the averaged current when a gate potential is 
used to move the levels of the dot. These peaks were 
observed experimentally^ and are associated with ab- 
sorbtion and emission proces of energy quanta from the 
pumping fields. Theoretically they were obtained first in 
the Master equation framework^ In the Keldysh formal- 
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ism a calculation of the stationary current in the WBL 
approximation was also presentent in Ref. l28l Here we 
want to check wether the satellite peaks appear also in 
the first first pumping cycles, namely in the transient 
regime. To this end we consider the dependence of the 
pumped current for a two site turnstile on a gate po- 
tential V g (see the Hamiltonian given in the introduc- 
tion) which shifts the two levels of the isolated system. 
Since we are interested in observing satellite peaks asso- 
ciated to one level only we take the parameters such that 
Vi r < u> < S2, where 6 = Ei~ E2 = 1 is the level spacing. 
Fig. 8a shows the averaged current associated to the 2nd 
pumping cycle One notices as once two satellite peaks 
located on each side of main resonant peaks (at V g = ±1 
the levels of the isolated dot are shifted to zero, threfore 
they are located in the middle of the bias window). Also, 
the distance between the sattelite peaks and the associ- 
ated main peak equals roughly the frequency lu = 0.65 
and confirms the absorbtion/emission picture. There are 
however several aspects due to the transient regime in 
which these peaks appear, i) As shown in Fig. 8b the av- 
eraged current on the first pumping cycle does not display 
satellite peaks and is mostly negative, ii) Also, in Fig. 8a 
the first sattelite peak has a negative value. This can be 
understood by looking at Fig 8c which gives the 3D plot 
of the pumped current as a function of time and gate 
potential. Clearly, around V g = —1.65 there is a max- 
imum positive current at the beginning of the pumping 
cycle. In this range the system effectively pumps charge 
to the right lead via photon assited tunneling involving 
the highest level of the dot. This process is rather weak 
because in the transient regime the occupation of the lev- 
els below the bias window is not complete. This is why 
after a short pumping regime the system absorbs charge 
from the lead and threfore J r (t) < 0. As the first level en- 
ters the bias window a main peak appears around V g = 0. 
Similar description can be made for the contribution of 
the second level, iii) For V g > 0.3 the average current 
is positive; the photon-assisted tunneling contribution to 
the pumping process is amplified and the current van- 
ishes as the lowest level of the dot is pushed upwards. 

Our calculations show that photon assisted tunneling 
takes place also in the transient regime. 



IV. CONCLUSIONS 

We have studied the turnstile pump regime of 
a few level noninteracting quantum dot within the 
non-equilibrium Green-Keldysh formalism for time- 
dependent transport. The out-of-phase oscillating bar- 
riers coupling the system to the leads are described by 
time-dependent hopping terms. In the numerical calcu- 
lations we have considered a train of trapezoidal pulses 
that mimick the configuration used in the experimental 
work of Kouwenhoven et al. 6 . Our approach includes 
explicitely the starting time of the pumping cycles and 
therefore captures the transient behavior of the time- 
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FIG. 8: (Color online) (a) The averaged current J r ,2 as a func- 
tion of the gate potential for a two sites turnstile shows satel- 
lite peaks associated to photon-assisted tunneling during the 
second pumping cycle, (b) The averaged current correspond- 
ing to the first pumping cycle J Tt i does not contain satellite 
peaks, (c) The current J r as a function of time (only the 
second cycle is shown) and gate potential. The discussion is 
made in the text. Other parameters lu — 0.65, vi = v r = 0.35, 
fit = 0.2, fi r = -0.2, kT = 0.0001. 



dependent and period-averaged current. To our best 
knowledge, such a calculation is presented here for the 
first time. We identify basically two stages of transport. 
The system experiences first a transient regime with poor 
pumping efficiency. This is due to the fact that the dot 
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rather absorbs charge from the leads in order to populate 
the levels located below the bias window. The number 
of pumping cycles in which the transient features should 
be observed in future experiments with many level turn- 
stiles increases if the bias window contains only the high- 
est levels. In the second stage the occupation number 
of the dot and the currents oscillate with the pumping 
period. We show that at low frequency and strong cou- 
pling to the leads an integer or half-integer number of 
electrons are pumped, depending on the number and the 
location of the levels within the bias window. In the high- 
frequency case the pumped charge is rather small even 
if additional contributions appear due to scattering pro- 
cesses involving energy sidebands. We show that satellite 
peaks due to the photon-assisted tunneling appear also 



in the transient regime. The present analysis could be 
extended to two-dimensional systems in order to discuss 
magnetic field effects. Also, different types of potentials 
(e.g. harmonic or damped pulses) can be considered. 
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